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Abstract —In this paper, a rather compact differential 
equation governing the bending behavior of a 
magneto-electro-elastic (MEE) rectangular thin plate is 
introduced, in particular, the exact solutions for the deformation 
response of laminated BaTi03-CoFe204 composites subjected 
to certain types of surface loads are analytically obtained. Due to 
the omission of the transverse shear deformation and rotatory 
inertia assumed in Kirchhoff thin-plate theory, the governing 
equation can accordingly be expressed in terms of the transverse 
displacement only. As a result, the structural characteristics 
such as elastic displacements, electric potential and magnetic 
induction for a magneto-electro-elastic (MEE) rectangular plate 
can be carried out in a theoretical approach. For a laminate 
MEE composite, the material constants can be uniquely 
determined by the volume-fraction (v.f.) of the piezoelectric 
constituent BaTi03, and are tabulated with 25% offset of the 
volume-fraction. According to the specified boundary conditions 
imposed on the MEE thin plate, the deformation variations with 
closed-circuit electric restriction are evaluated analytically in 
the present study. The results obtained in this paper by using the 
proposed model can be shown to have good agreements with the 
other available research works, however, with the advantage 
that the present study indeed provides a much simpler way in 
seeking the analytic solutions for the interactively coupled 
quantities of a layered MEE medium. 

Index Terms —Layered structures, Mechanical properties, 
Analytical modelling 


I. Introduction 

It is well-known that structure made of piezoelectric 
materials can produce voltage when external stress is applied, 
nevertheless, it will also induce stress when voltage difference 
is implemented across the structure. Parallel to the 
piezoelectricity, piezomagnetizm can intrinsically be 
characterized by interactive couplings between system's 
magnetic polarization and mechanical strain, i.e., in a 
piezomagnetic medium, one may induce a spontaneous 
magnetic moment by applying physical stress, or reversely 
obtain a physical deformation by applying a magnetic field. 
Composites structures consist of layered piezoelectric and 
piezomagnetic components, which can be referred as one kind 
of the magneto-electro-elastic (MEE) material, possesses the 
ability to convert energy of magnetism, electricity or elasticity 
into another form. The MEE material also exhibits a specific 
magneto-electric effect which is not appeared in a 
single-phase piezoelectric or piezomagnetic material. 
Surprisingly, in some cases the magneto-electric effect of 
MEE composites can even be obtained two orders larger than 
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that of a single-phase magneto-electric material with highest 
magneto-electro coefficient [1]. Due to these multiphase 
properties, the MEE material have been found increasing 
applications in making efficient smart and intelligent 
structures such as magnetic field probes, electric packing, 
acoustic, hydrophones, medical ultrasonic imaging and so on. 

For the past few years, systematic investigations either in 
determining the material coefficient of such new type material 
or in analyzing the static or dynamic behavior under certain 
external conditions are vigorously proposed by many 
professional engineers and scientists. Li [2] studied the 
average magneto-electro-elastic fields in the multi-inclusion 
embedded in an infinite matrix, and estimated the effective 
magneto-electro-elastic moduli of piezoelectric- 
piezomagnetic composites for both BaTi03 fiber reinforced 
CoFe204 and BaTi03-CoFe204 laminate. It has been shown 
that the magneto-electro coupling demonstrated by 
magneto-electro coefficients vary with the volume faction of 
BaTi03 and have opposite signs for fibrous and laminated 
composites, respectively. Pan [3] obtained the exact solution 
for three dimensional, anisotropic, linear 
magneto-electro-elastic, simply-supported and multilayered 
rectangular plates under static loadings. It is stated that even 
for relatively thin plate, responses from an internal load are 
quite different to those from a surface load. Chen et. al. [4] 
establish a micro-mechanical model to evaluate the effective 
properties of layered magneto-electro-elastic composites and 
the linear coupling effect between elasticity, electricity and 
magnetism of the MEE composite is accordingly derived. In 
their study, numerical results for a BaTi03-CoFe204 
composite with 2-2connectitivi are obtained, and the 
dependences of magneto-electro-elastic coefficients, the 
so-called product properties, of the composite on the volume 
fraction of BaTi03 are clearly depicted. Wang and Shen [5] 
extended their previous works on piezoelectric media to study 
the general solution of three-dimensional problems in 
transversely isotropic magneto-electro-elastic media through 
five newly introduced potential functions. Chen and Lee [6] 
simplified the governing equations of the linear theory for the 
magneto-electro-thermal-elastic plate with transverse 
isotropy by introducing two displacement functions and stress 
functions. In their study, two new state space equations are 
established while selecting certain physical quantities as the 
basic unknowns. Wang et. al. [7] performed a state vector 
formulation for the three dimensional, orthotropic and 
linearly magneto-electro-elastic multiple layered plate and 
expressed the basic unknowns by collecting not only the 
displacement, electric potential and magnetic potential but 
also some of the stresses, electric displacements, and 
magnetic induction. A boundary integral formulation for the 
plane problem of magneto-electro-elastic media are 
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performed by Ding and Jiang [8] using strict differential 
operator theory. They obtained the fundamental solutions for 
an infinite MEE plane in terms of four harmonic functions 
which satisfied a set of reduced second order partial 
differential equation for distinct eigenvalues case. Guan and 
He [9] derived the fundamental equation for the plane 
problem of a transversely isotropic magneto-electro-elastic 
media by applying the Almansi’s theorem and expressed all 
physical quantities by four harmonic functions for distinct and 
non-distinct cases. 

In this study, a rather simple analytic solution for the 
deformations of the magneto-electro-elastic (MEE) 
rectangular plate under certain type of applied loads acting on 
the top surfaces are derived. By imposing the Kirchhoff 
thin-plate hypothesis on the plate constituent, the governing 
equation in terms of only the transverse displacement of the 
plate can be obtained and therefore a rather compact form 
indicating the multiple effects between elasticity, electricity 
and magnetism of the plate can be successfully presented. The 
MEE plate is chose to be made of the two-layered 
BaTi03-CoFe204 laminate, which can be thought as a 
transversely isotropic magneto-electro-elastic medium and 
the material coefficients for such continuum can be expressed 
uniquely by introducing the volume-fraction of BaTi03 in the 
layered composite. The corresponding deformation analysis 
regarding the elastic displacements, electric potential and 
magnetic induction of the MEE thin plate is evaluated through 
the formulation mentioned in this study. Some comparisons 
with previous literatures are made and great agreements are 
reached which directly validate the proposed simplification 
for the MEE modeling. 


II. Formulations 


For a Magneto-Electro-Elastic (MEE) thin plate made of 
two-layered BaTi 0 3 -CoFe 2 04 laminate, the fundamental 
Kirchhoff hypothesis for the small-deflection of simple 
bending problem can be applied, and thus the transverse shear 
deformations and rotary inertias can be neglected, also the 
transverse shear strains are negligible. In accordance with the 
assumptions that in-plane electric fields and magnetic fields in 
a very thin medium can be ignored [10], that is, only the 
transverse electric, E 3 , and magnetic field, H 3 , are under 

consideration, the following governing equation for the 
bending problem can be found [11] 

|dV 4 h'+ £V 4 h'+MV 4 w1 = P(x,y ), (1) 


where D = 


Ci i h 


Ml 


12 


represents the plate rigidity, 


e 3] h A, 

E = — - - denotes the effective rigidities induced by 

12 A 

q 3] h 2 A ? 

electricity, M = —— -is the corresponding effective 

rigidities cause by magnetism and P(x , y) denotes the 


applied load acting on the top surface. Herein h is the plate 
thickness, A = £ 33 p 33 ~d^ » ^1 - “^ 33 ^ 33 ) an ^ 

A 2 = ~{s 33 q 3l ~d 33 e 3l ) are extra-defined parameters aiming 
to simplify the derivations; c t - , e- , q- , d- and ju- are the 

elastic, dielectric, piezoelectric, piezomagnetic, 

magnetoelectric, and magnetic constants, respectively. It 


should be noted the relation c n = c l2 +2 c 66 for transversely 

isotropic material is adopted in the above equation. As it can 
be also learned from Ref. [11], the reduced extended traction 
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Meanwhile, the following expression for the electric potential 
and magnetic potential can be derived, 


~ = -■^-zV 2 w + ^ 1 (x, y), (5) 

dz A 

^- = -^yzV 2 w + y/\(x,y), (6) 

dz A 

where q\(x,y) and y/ 1 (v, y) represent the variations of 

electric field and magnetic field in the thickness direction 
while the plate is under deformation, and are both 
independent of z variable. 


In seeking for the solution to Eq. (1), we can assume the 
following expression for the transverse deflection of the MEE 
plate, 

oo oo 

w(x, y ) = 'E J EjA mn X m (x)Y n (y) , (7) 

m =1 77=1 


where X m (x) and Y n (y) are the homogeneous solutions of 

Eq. (1) and can be determined according to the assigned 
boundary conditions. Some mode shapes and corresponding 
eigenvalues with respect to commonly seen boundary 
conditions are tabulated in Table 1 as a reference. It should be 
noted that the mode shapes X m (x) and Y n (y) not only 


satisfy the corresponding boundary conditions but also 
possess the intrinsic orthogonality, i.e., we will have 


Jo 0, if m ^ M 
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( 8 ) 


and 

[ Y n (y)Y N (y)dy 

J 0 


Y n (y )|, if n = N 
0, if n*N 


Y„(y ) S, 


nN ? 



where L x and L v denotes the plate lengths along x and y 

directions. After X m (x) and Y n (y) are determined, we can 

further expand the applied load on the top surface of the plate 

into the generalized double Fourier series as 

00 00 

p O,y) = 'EXj p mn x m (x)Y n (y), (10) 

777=1 77=1 

in which the Fourier coefficient p mn can be determined as 
follows 
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P mn 


1 rU rL 

XJx) K(y)||Jo Jo 


J J P(x, y)X m (x)Y n (y)dxdy . ( 11 ) 


By substituting Eq. (7) and Eq. (10) into Eq. (1) and 
collecting the constant term, we can have the following 
equation 

(D + E + M) ( a A m + 2+ Pi j A mn = p mn , (12) 

where a m and p n are the corresponding eigenvalues 

associated to the specific boundary conditions. Furthermore, 
the magnitude of the transverse deflection, A mn , can be 

determined, i.e.. 

Am 


Case I: Closed-circuit, i.e., (f)(x, y , ±h /2) = 0 and 

y/(x , y,±h/2) = 0 

By carrying out the anti-derivatives for Eq. (5) and Eq. ( 6 ), we 
can have the following expressions for electric potential and 
magnetic potential as 

<p(x, y,z) = ~7 z 2 V 2 w(x, y) + zq\ (x, y) + %(x,y) (17) 


2A 
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y/(x, y , z) = ~ -f z 2 V 2 w(x, y) + zy/\ (x, y) + y/ 0 (x, y) (18) 

2A 

since y , ±h / 2 ) = 0 , after solving the above two 
equations, one can get 
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due to applied load 
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and the exact solution for the mechanical displacements along 
x- and y- directions can be expressed as 


W (*,;y) = 0 and %(x,y) = ^X(h/2) V w(*,y), 

thus the exact solution for the electric potential due to applied 
load can be expressed as 

<p(x, y,z) = -2 L z 2 V 2 w(x, y) + A- (h / 2 ) 2 V 2 w(x, >>). (19) 

2A 2A 

By the same token, we can have the following exact solution 
for the magnetic potential due to applied load as 
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Table 1: Mode shapes and the corresponding eigenvalues for 
specified boundary conditions 


Case II: Open-circuit, i.e., (x, y , ±h / 2) = 0 and 

B z (x, y,±h/2) = 0 

Substituting Eqs. (5)-(6) and Eqs. (15)-(16) into Eqs. (3) and 
(4), the solutions for the electric displacement and magnetic 
induction in thickness direction can be expressed as 
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As for the electric boundary conditions and magneto 
boundary conditions specified on the top and bottom surfaces 
of the plate, there are normally two kinds of cases to be 
discussed. 


B z (x,y,-hl 2 ) = - -<731+^33 y + ^33 yy (^/ 2 )V : 2 w-d 33 (p x (x,y)-ju 33 \f/ l (x,y) = 0 . (26) 

V A A 

In viewing of Eq. (23) through Eq. (26), it can be found that 
these four equations cannot provide any information about 
determining the expressions for cp x (x, y) and y/ x (x, y), nor 

for the terms (p 0 (x,y) and y/ 0 (x,y) . This is because the 

variations of electric displacement and magnetic 
displacement along thickness direction are vanished under the 

DD dB 

proposed methodology, i.e., —— = —— = 0 , therefore it is 

dz dz 

expected to see the linear dependency of these two quantities 
on the z variable. In that sense, open-circuit boundary will 
result in a trivial solution with the displacement components 
approaching to zero and makes the expressions of cp x {x, y) 

and y/ x (x,y) being negligible as can be detected from Eqs. 
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(21) and (22). As a result, two more boundary conditions will 
be needed in order to resolve the expressions for the terms 
<Po (x, y ) and y/ 0 (x, y) . For that reason, in this study, we 

don’t specifically focus on the open-circuit case and only 
adopt the closed-circuit boundary condition as our numerical 
examples discussed in the next section. 

It should be noted that the in-plane electric fields and 
magnetic fields can be ignored if the plate thickness is very 
small (e.g. h <10L x ), and only the transverse electric field, 


E 3 , and magnetic field, H 3 , are related to the electric 
potential cp and magnetic potential y/ in the follow form 
according to the Maxwell’s equations. 

E =-^- = ^-zV 2 w-(p i {x,y), (27) 

dz A 

H 7 =-^- = ^yz^ 2 w-y/ x (x,y), (28) 

dz A 

in which the terms cp x (x, y) and y/ x (x, y) can be determined 

by considering the closed-circuit case on the plate surfaces. 
By substituting Eqs. (5)-(7) into Eq. (2), the stress 
distributions are respectively 

a x (x,y,z) = -Cn^zK^Y^-C^zK^Yfy) 
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'31 A + ^31 A 

A A 
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+e 3i <p l (x,y) + q il y/ i (x, y) 
a y (X, y, z ) = -CnAnn^m (x) Y „(y) - Cl 1A m 


1 Ay 

31 A + ^ 31 A 

A A j 


f A 
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A mn z{x:(x)Y n (y) + X m (x)Y;(y)), (30) 


+e 3l <p l {x,y) + q il i// l (x,y) 

T xy- = - 2c 6eA nnZX'nWY^iy) . (31) 

The above equations provide the analytic solutions to the 
physical quantities of a magneto-electro-elastic thin plate 
subjected to surface applied load and can be decided 
according to the transverse deflection w(x, y) and potential 


parameters <p(x, y, z) and y/(x, y, z ). 


III. Numerical Demonstrations 

In this section, some examples based on the proposed 
model for the thin magneto-electro-elastic rectangular plate 
subjected to external loads and closed-circuit electric 
restrictions are presented. For the first a few examples, the 
author demonstrates some special cases which have been 
verified in previous research works in order to validate the 
present study and tell the differences with other available 
literatures. 

Mode shape validation. Since the mode shape 
corresponding to the free vibration of the MEE plate is of 
interest due to its importance on the nature and extent of the 
electro-magneto-mechanical coupling, we first check the 
validation of fundamental mode shapes for a piezoelectric 
square plate made of commonly used PZT4 material by 
utilizing the present model. In order to compare the results 
with those proposed by Heyliger & Saravanos [13], the 
dimensions of the plate are chosen to be L x / h - 10 , L x - L y 

and h = 0.01 m , and the boundary conditions are assumed to 
be simply-supported around four edges with closed-circuit on 


the top and bottom surfaces, i.e. 
^(x, y, -EL / 2) = y,-EL / 2) = 0 . The thickness is 

normalized as Z = z/ H , the displacements and electric 
potential are normalized with respects to their own largest 
values along thickness direction, i.e., TT = w/max(w) , 

U=u/max(u ) and O = 0/ma x(u) . The normalized 

distributions for both the transverse and in-plane mode shapes 
along with the electric potential for a PZT4 square thin plate 
are depicted in Fig. 1, as we can detect from this figure, it is 
almost identical to Fig. 2(a) in the paper conducted by 
Heyliger & Saravanos [13]. Therefore, good agreement has 
been reached. 



Fig 1. Thickness distributions for a PZT4 thin plate with 
dimension L x l h -10 , L x = L and h = 0.01 m subjected to 

simply-supported and closed-circuit boundary conditions. 

Example 1. The second example is an electro-elastic single 
layer rectangular plate made of purely piezoelectric BaTi0 3 . 
The length-to-thickness ratio and the width-to-thickness ratio 
are chosen to be L v / H - 10 and L v / H = 5 in order to 

match the dimension setting in the paper presented by Chen 
and Lee [6] although this doesn’t assure the thin-plate 
requirement to be satisfied. All the surfaces of the plate are 
assumed to be traction free except the bottom surface, on 
which a sinusoidal loading 

A P - <T 0 sin (nnx / L x ) sin(m^y / L y ) with amplitude 

cr 0 =lN/m and mode number m = n = 1 is applied. In 

particular, the deformation responses of the plate are 
calculated at the fixed horizontal coordinates 
(x 0 , y 0 ) = (0.75L r ,0.25L v ) . Fig. 2(a)-(f) indicates the 

variations along thickness direction of the elastic 
displacements w , u x , u y , electric potential ^(x, y, z) , 

magnetic potential y/(x, y, z) , electric displacement D z and 

magnetic displacement B , for such a plate caused by a 

sinusoidal loading on the top surface with all edges 
simply-supported. As we can recognize from these figures, 
the shear deformation is linearly dependent on the transverse 
deformation and the electric potential reveals quadratic 
variation along z coordinate. These two phenomena can be 
also found in both the papers studied by Pan [3] and Chen [6], 
see Fig. 1 and Fig. 2 in their studies respectively. It should be 
noted that the electric potential and magnetic potential in this 

DD cB 

case are identically zero due to the fact that —- = —- = 0 

dz dz 
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and close-circuit electric restriction is adopted in this paper 
instead of open-circuit one. 


(a) transverse deformation w(x,y) 


(b) shear deformations 
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Fig. 2. Variations along thickness direction of the elastic 
displacements, electric potential, magnetic potential, electric 
displacement and magnetic displacement for a single layer 
BaTi03 plate caused by sinusoidal loading on the top surface 
with all edges simply-supported. 


Table 2. Material constants for both laminate and fibrous 
MEE multiphase composites, partially cited from Buchanan 
[ 12 ] 


Material 

constants 

PZT-4 

CoFe 2 0 4 

Bi-layered BaTi0 3 -CoFe 2 0 4 laminate 

BaTi0 3 

v.f. 

N.A. 

0% 

25% 

50% 

75% 

100% 

<11 

138.5 

286 

250 

215 

188 

166 

Q: 

77.3 

173 

141 

112 

90 

77 


73.6 

170 

139 

110 

88 

78 

% 

114.7 

269.5 

245 

225 

205 

162 

Q 4 

25.6 

45.3 

45 

45 

45 

43 

*31 

-5.2 

0.0 

0.0 

0.0 

0.0 

-4.4 

*33 

15.08 

0.0 

0.1 

0.25 

0.4 

18.6 

*15 

12.72 

0.0 

3.0 

6.0 

8.8 

11.6 

*11 

1.306 

0.08 

3.4 

6.3 

8.8 

11.2 

*33 

1.115 

0.093 

0.2 

0.25 

0.3 

12.6 

Mi 

5.0 

-590 

-4.45 

-2.95 

-1.55 

5 

/'33 

10 

157 

0.35 

0.20 

0.05 

10 

<h\ 

0 

580.3 

80 

10 

5 

0.0 

?33 

0.0 

699.7 

100 

30 

15 

0.0 

<hs 

0.0 

550 

410 

255 

130 

0.0 

Mi 

0.0 

0.0 

-27500 

-35500 

-27000 

0.0 

M3 

0.0 

0.0 

-24.5 

-24.8 

-23.5 

0.0 


(Cy in 10 9 N/m 2 , e y in C/m 2 , q v in N/Am 2 , in 10’ 9 C 2 /Nm 2 ,and in 10^Ns 2 /C 2 


and Wj, in lO'^Ns/VC) 

Example 2. In this example, a two-layered 
magneto-electro-elastic plate made of equally-placed 
BaTi 0 3 -CoFe 2 04 laminate is then presented. The 
piezoelectric BaTi0 3 is placed on the bottom layer 


(-H / 2 < z < 0 ) whereas the magnetostrictive CoFe 2 0 4 is on 
the top layer (0< z < H /2), the material constants for this 
kind of MEE plate can be found in Table 2 with the 
volume-fraction set to be 50% for the lamination case. The 
dimensions are L v = 1000 m , L = 20 m , h = 1 m , 

S = L / h - 20, n = 1 and the related physical quantities are 
normalized as follows with all of them still remain 
dimensional, u/S 3 , v/S 3 , wxlOO/S 4 , <J XX /S 2 , cr yy /S 2 , 

<7 xy /5 2 , (,MS 2 , If/IS 2 , DJS , and BJS . The 

deformation variations of the magneto-electro-elastic plate 
due to external load, A P(x, y) = P 0 sin ny / L y , applied on the 

top surface with magnitude P 0 =lN/m are presented in Fig. 
3 at the location (v 0 , y 0 ) = (L x / 2,574). Fig. 3(a)-(h) are the 
variations of the elastic displacements components w(x,y) 
and u y (x 0 ,y 0 ,z ) , electric potential (f>(x 0 ,y 0 ,z) , magnetic 

potential y/(x 0 ,y 0 ,z) , electric displacement D„ , magnetic 

displacement ZT , the normal stress components 

cj n (x 0 ,y 0 ,z) and a 22 (x 0 , y 0 , z) along the thickness 

direction with boundary condition to be simply-supported 
around four edges. From these figures, we can obviously 
observe the interactive behavior between the piezoelectric, 
piezomagnetic and magnetoelectric effects for a 
magneto-electro-elastic (MEE) plate under mechanical 
applied load. If the MEE plate is assumed to be very thin, the 
deformation variations for both the electric and magnetic 
potentials reveal a quadratic dependence on the thickness 
variable, however, for the other related quantities such as 
shear deformations, electric and magnetic displacements as 
well as the stress distributions, linear dependence can be still 
detected. 

From now on, the bi-layered BaTi0 3 -CoFe 2 0 4 laminate 
working as a magneto-electro-elastic plate will be discussed 
based on the Kirchoff hypothesis and the thin-plate theory. 
The simplified governing equation, Equation (1), will be 
examined subjected to various kinds of surface applied load, 
and the deformation behavior with different boundary 
conditions imposed on the plate will also be inspected. Due to 
the simplicity of the proposed model, the higher mode 
response can be easily performed; therefore, the deformation 
variation of the MEE plate with respect to different mode 
numbers will be carried out as a reference. 

The dimensions of the MEE plate is set to be 
L x xL xH = lm*lm*0.05m unless otherwise mentioned, 

however, any kind of dimensions can be applied as long as the 
span-to-thickness ratio is satisfying the requirement for 
thin-plate theory, i.e., L x 7 H > 10 . The plate surfaces are 

assumed to be traction free except on the top or bottom 
surface, on which a z-direction surface load is applied. The 
external load can be of any type possibly occurs in the study 
of MEE plate, however, in order to observe the variation of 
the deformations, three commonly seen static forces are 
performed in this paper. They are 

(1) uniform load, i.e., P(x,y) — P 0 , 

(2) distributed load, i.e., P(x,y) = P 0 X M (x)Y N (y ) and 

(3) concentrated load, i.e., P(x,y) = P 0 S(x-x 0 ,y - y 0 ). 
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(a) transverse deformation w(x,y) shear deformation v(x Q ,y 0 ) 



(c) electric potentail <Kx 0 ,y 0 ) (d) magnetic potentail y(xO,yO) 




(e) electric displacement D 2 (x QI y 0 ) 


(f) magnetic dispalcement B z (x 0 ,y 0 ) 


MEE plate with simply-supported around edges and surface 
distributed load AP = sin(7rx/L x )sin(7ry/L y ) with respect 

to various volume fraction of BaTi03. As it can be observed, 
the concavity for the pure piezoelectric BaTiC^ plate 
(v.f.=100%) is uniquely different from the other MEE 
laminates with a much higher magnitude. Although the 
dependence with respect to z variable remains quadratic, the 
electric potential for the MEE lamination with volume 
fraction other than 100% reveals a negative sign along the 
thickness direction instead of a positive one in the pure 
BaTi 03 case. However, the magnitudes placed in order are 
25%, 75% and 50%, 0% for the pure piezomagnetic CoFe204 
case stands zero with no piezoelectric effect due to the applied 
load, which is quite reasonable. As for the magnetic 
potentials, quadratic relation can also be verified and the 
concavity for each case seems to be consistent. The maximum 
magnitude occurs in the 25% laminate followed by the 75%, 
50% and 0% cases, 100% for pure piezoelectric BaTi 03 case 
stays nil with piezomagnetic effect being vanished. 



4(x 0 ,y 0 ,z) 


Fig. 3. Variations of the elastic displacements, electric 
potential, magnetic potential, electric displacement, magnetic 
displacement and normal stress components along the 
thickness direction for a two-layered magneto-electro-elastic 
BaTio3-CoFe204 plate caused by external loading on the top 
surface with all edges simply-supported. 

The geometric boundary conditions imposed on the MEE 
plate can be any one of the combinations for free, clamped 
and simply-supported edges; however, in this study, only the 
boundary conditions of simply-supported around, clamped 
around and cantilever plates are conducted as numerical 
examples. It should be noted that the deformation responses 
for the MEE plate under static loads will be calculated at a 
fixed horizontal coordinate (x 0 , y 0 ) , however, due to the 

different plate characteristic with respect to different 
boundary conditions, the horizontal coordinates will be 
various according to the corresponding boundary conditions. 
That is, for SSSS plate location is chosen to be at 
( Xo ,y 0 ) = (0.5L x ,0.5L v ) , for CCCC plate at 

(x 0 ,y 0 ) = (0.5 L y , 0.5L } ,) and for CFFF plate at 

(x 0 ,y 0 ) = (L x ,0.5L y ). 

Fig. 4 (a)-(d) are the deformation variation of electric 
potential, magnetic potential, electric displacement and 
magnetic displacement for the laminated BaTi 0 3 -CoFe 2 04 



v<W> 



Fig. 4-1. Variation of (a) electric potential (b) magnetic 
potential for the lamination MEE plate with simply-supported 
BCs under surface distributed load. 
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D z <vyo’ z) 



B z <w> 



xIO 23 


Fig. 4-2. Variation of (c) electric displacement (d) magnetic 
displacement for the lamination MEE plate with 
simply-supported BCs under surface distributed load. 

For the MEE laminate presented in Fig. 4, the electric and 
magnetic displacements, also called fluxes, can be found to be 
linearly dependent on the z variable and the slope varies 
slightly with respect to the volume fraction of the MEE plate. 
However, if we near watch the magnitudes, we can find that 
all of them are pretty slim. This is due to the nature of the 
material parameter itself, also because of the close-circuit 
restriction we impose on the model. Nevertheless, if the 
electric and magnetic boundary conditions are chosen 
otherwise, the magnitudes for electric and magnetic 
displacements may be changed in a different way. 
Whatsoever, in the present study, only the close-circuit MEE 
thin plate is considered, therefore the author hereafter leave 
out all the presentations and discussions for the electric 
displacement and the magnetic displacement due to their 
insignificant effects. 

Since the deformation variations for the laminate MEE 
plates with clamped-around edges and under uniform applied 

load A P = 1 N/m are quite similar to those behaviors 

presented in Fig. 4, the author therefore skip these two case 
and directly go to the cases for the cantilever laminate and 
fibrous MEE plate. Also the deformation variations for the 
electric and magnetic displacements are neglected due to their 
small magnitudes both approaching to zero. 

Fig. 5 (a)-(b) are the deformation variation of electric 
potential, magnetic potential versus volume fraction for a 
cantilever laminated BaTi 0 3 -CoFe 2 04 plate subjected to unit 


ISSN: 2394-3661, Volume-3, Issue-4, April 2016 

impulse force acting on the location (x 0 ,y 0 ) , i.e., 

A P = 1* S(x- x 0 ,y- y 0 ) • The deformation behavior of a 

cantilever MEE plate is quite similar to the corresponding one 
of a simply-supported MEE plate except for the concavity. 
Owing to the different characteristics of MEE plate with 
different boundary conditions and subjected to different types 
of applied load, the sign change on the concavity is reasonable 
and expected. 


4(x 0 ,y 0 ,z) 






Fig. 5. (a)-(b) Deformation variation of electric potential and 
magnetic potential versus volume fraction for a cantilever 
laminated BaTi 0 3 -CoFe 2 04 square plate subjected to unit 
impulse force acting on the location (x 0 ,y 0 ) , i.e., 

A P = l-S(x-x 0 ,y-y 0 ) . 

The deformation variation of electric potential and 
magnetic potential for the 50% laminated BaTi0 3 -CoFe 2 0 4 
MEE plates subjected to various loads with all edges 
simply-supported are demonstrated in Fig. 6 (a)-(b). And the 
corresponding deformation variations for the laminate MEE 
composites with all edges clamped are depicted in Fig. 7 
(a)-(b). Followed by Fig. 8 (a)-(b), the same illustrations for 
the MEE cantilever plates are provided. As we can see from 
these figures, the concentrated applied load always stimulate 
much stronger deformation as we expected and followed by 
the distributed applied load, uniform applied load seems to be 
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Fig. 6. (a)-(b) Deformation variation of electric potential and 
magnetic potential for the 50% laminated BaTi0 3 -CoFe 2 0 4 
MEE plate subjected to various loads with all edges 
simply-supported. 
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Fig. 7. (a)-(b) Deformation variation of electric potential and 
magnetic potential for the 50% laminated BaTi0 3 -CoFe 2 0 4 
MEE plate subjected to various loads with all edges 
simply-supported. 
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Fig. 8. (a)-(b) Deformation variation of electric potential and 
magnetic potential for the 50% laminated BaTi0 3 -CoFe 2 0 4 
MEE cantilever plate subjected to various loads. 

Finally, in order to see the influence of mode orders, the 
deformation variation of electric potential and magnetic 
potential versus mode number for the MEE laminate under 
certain applied load are given in Fig. 9 (a)-(b) with all edges 
simply-supported, in Fig. 10 (a)-(b) with all edges clamped, 
and in Fig. 11 (a)-(b) for the cantilever plate. It should be 
noted that some of the mode orders contributes nothing on the 
plate deformations because of maybe the symmetry of the 
mode shape or the location we calculate at. 


(a) (b) V (x 0 ,y 0 ) 




Fig. 9. (a)-(b) Deformation variation of electric potential and 
magnetic potential versus mode number for the 50% 
laminated BaTi0 3 -CoFe 2 0 4 MEE simply-supported plate 
subjected to sinusoidal loading. 




(b) v(^.y 0 ) 
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Fig. 10. (a)-(b) Deformation variation of electric potential 
and magnetic potential versus mode number for the 50% 
laminated BaTi0 3 -CoFe 2 0 4 MEE clamped plate subjected to 
sinusoidal loading. 
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*<Vo 



(b> v(vy 0 > 



Fig. 11. (a)-(b) Deformation variation of electric potential 
and magnetic potential versus mode number for the 50% 
laminated BaTi0 3 -CoFe 2 0 4 MEE cantilever plate subjected 
to sinusoidal loading. 


IV. Conclusion 

The closed form solutions for the bending problem of a 
bi-layered BaTi 0 3 -CoFe 2 04 composite are derived based on a 
new invented governing equation for magneto-electro-elastic 
(MEE) rectangular thin plate, in particular, the elastic 
displacements, electric potential and magnetic induction for a 
magneto-electro-elastic (MEE) laminate are implemented 
analytically. 

It has been shown that the material coefficients for the 
MEE constituent vary a lot according to the volume-fraction 
of BaTi0 3 it contains. The deformation variations for the 
MEE thin plate with closed-circuit electric restriction are 
evaluated with respect to various boundary conditions, and 
the effects of the volume-fractions are investigated in detail. It 
can be found that the shear deformation is linearly dependent 
on the transverse deformation, whereas the electric and 
magnetic potential are both of quadratic variation along the 
thickness direction. In addition, the deformation behavior for 
a single phase material can be found to be quite different from 
the multiphase one in either the magnitude or the sign it is 
induced by. 

The present study provides some commonly seen examples 
for the magneto-electro-elastic (MEE) rectangular thin plate 
under the action of 3 kinds applied loads, and offers the 
discrepancy on the deformation variation of electric potential 
and magnetic induction with respect to various typical 
boundary conditions. This work proposed a much easier and 
systematic way to seek for the analytic solutions for the 
deformation characteristics of a bi-layered MEE thin plate, 
and should be of interest to someone devoted on the practice 
of structure design with the fully coupled medium. 
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